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Abstract 



c/2 ' We study stable axially and spherically symmetric spatial solitons in plasma 

1—1 ■ with diatomic ions. The stability of a soliton against the collapse is provided 

. . by the interaction of induced electric dipole moments of ions with rapidly os- 

Q ! dilating electric field of a plasmoid. We derive the new cubic-quintic nonlin- 

' c/3 \ ear Schrodinger equation which governs the soliton dynamics and numerically 

^ 1 solve it. The stability of solitons is analyzed using the Vakhitov-Kolokolov 

P^' criterion. Then we discuss the possibility of implementation of such plas- 

moids in a realistic atmospheric plasma. In particular, we suggest that ax- 
ially and spherically symmetric Langmuir solitons is a theoretical model of 
QQ \ long-lived atmospheric plasma structures. 

^ ■ Keywords: spatial Langmuir soliton, induced electric dipole moment 

O ■ long-lived atmospheric plasma structure 

cn 

^ 1. Introduction 

_^ Stable spatial solitons are observed in the studies of optical phenom- 

^ \ ena jl|, in sohd states physics |2|, and in the plasma research |3|. Typically 

the stability of a soliton is provided by a certain nonlinearity. For example, 
spatial plasma solitons, described in frames of the classical electrodynamics, 
can exist due to the combined action of electron-ion and electron-electron 
nonlinear interactions. The former one was found in Ref. j^ to be focus- 
ing, whereas the latter interaction can be defocusing |5|, |6|. Recently it was 
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suggested [7| that, if one accounts for the additional quantum pressure of 
electron gas, it may explain the appearance of spatial Langmuir solitons in 
dense plasmas. 

In the present work we stall study the existence of stable Langmuir 
solitons in a plasma with ions possessing induced electric dipole moments 
(EDM). We shall demonstrate that the interaction of ion's EDM with a 
rapidly varying electric field of a plasma oscillation results in the arrest of 
the Langmuir wave collapse. Thus the existence of a spatial soliton becomes 
possible. 

In our analysis we shall choose a plasma with diatomic ions, which cor- 
responds to a realistic atmospheric plasma. Thus our results may be applied 
for the theoretical description of long-lived plasma structures observed in the 
atmosphere [8|. It is interesting to notice that the role of EDM of charged 
particles for the explanation of the stability of atmospheric plasmoids was 
also discussed previously in Ref. |9|. 

This work is organized as follows. In Sec. [2] we consider the general 
description of nonlinear waves in plasma. We introduce the new ponderomo- 
tive force, associated with EDM, which acts on the ion component of plasma. 
Then we derive a system of nonlinear equations for the electric field ampli- 
tude and the perturbation of the ion density. In Sec. [3] we reduce this system 
to a single nonlinear Schrodinger equation (NLSE) for the envelope of the 
electric field containing cubic and quintic terms. This equation is analyzed 
numerically for the case of a radial plasma oscillation. We show that both ax- 
ially and spherically symmetric stable solitons are possible. Then, in Sec. HI 
we discuss a possible application of our model of stable spatial solitons to 
the description of long-lived atmospheric plasmoids. Finally, in Sec. O we 
briefly summarize our results. The calculation of the permittivity of a gas 
with diatomic nonpolar molecules in an external electric field is presented 
in 



Appendix A 



2. Nonlinear Langmuir oscillations in plasma 

In this section we shall derive the basic nonlinear equations for the de- 
scription of Langmuir waves in plasma accounting for a ponderomotive force 
acting on nonpolar diatomic ions. 

If we study electrostatic plasma oscillations, i.e. when the magnetic field 
is zero, the motion of electron component of plasma obeys the following 



hydrodynamics equations: 

dn, 



' + V ■ KVe) = 0, 



dt 

^ + (v,VK = --E-^Vp, (1) 

at m Uem 

where He is the number density of electrons, Ve is the electron velocity, m 
is the mass of an electron, e > is the proton charge, E is the strength of 
the of the electric field, and p is the pressure. We should add to Eq. ([T]) 
the corresponding Maxwell and the Poisson equations for the electric field 
evolution, 

— - = 47re(neVe - rijVj), 

CM 

(V-E) = -47reK-ni), (2) 

where n, is the ion number density and Vj is the ion velocity. 

In the zeroth approximation only electrons participate in a plasma os- 
cillation, with the number density of ions being approximately constant 
m K. riQ = const. Thus we may present the electric field in the form, 

E = Eie"''^=* + E*e^'^'=* + ■ ■ ■ , (3) 

where Ue = \f^nehiQJm is the Langmuir frequency for electrons and Ei is 
the amplitude of an oscillation. 

In a realistic situation ions will also participate in a plasma oscillation. 
Thus the ion density becomes rij = no + n(r,t), where \n\ ^ tiq. It leads to 
the appearance of higher harmonics omitted in Eq. ([3]). The hydrodynamic 
equations for the description of the ions evolution has the form 



^ + V ■ (n^v,) = 0, 



dn 

'dt 

where M is the ion mass, (pes is the effective potential which obeys the 
following Poisson equation p: 

VVeff = 4vre (^^ + rlV'n) . (5) 



Here r^ = ^/Tj^nehi^ is the Debye length, T^. is the electron temperature, 
t/pm = |Eip/47r is the potential of the ponderomotive force which acts on a 
charged particle in a rapidly oscillating electric field ([3]). 

In Eq. (jl]) we include the volume density of the external force fpoi which 
acts on the ion component of plasma. Let us we suggest that ions are non- 
polar molecules which can acquire EDM pi = aijEj in an external electric 
field. Here {aij) is the polarizability tensor of an ion. If an ion is diatomic 
and possesses an axial symmetry, one can always reduce this tensor to the di- 
agonal form, {aij) = diag [a±, a±, a\\), where a± and on are transversal and 



longitudinal polarizabilities. Using the results of Ref. [10| and Eq. (lA.Sp . we 
get the ponderomotive force fpoi as 
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where e is the permittivity of the ion component of plasma, Tj is the ion 
temperature, (a) = (2a_|_ + a||)/3 is the mean polarizability of an ion, and 
Aa = a\\ — a^. 

Combining Eqs. ([I])-© we get the following nonlinear coupled equations 
for the electric field, 

iE + -a;erDV(VE) - ^nE = 0, (7) 

and for the perturbation of the ion density, 

d'^ 9^9\ 1 ^9,^,9 4 (Aa)^no„o,^,4 - , 

dt^ ' / 47rM ' ' 15 MT. ' ' ' ^ ^ 



where Cg = yT^/M is the sound velocity in plasma. Note that in Eqs. ([7]) 
and ([8]) we omit the index "1" in the amplitude of the electric field, i.e. 
El = E, in order not to encumber the formulas. 

It should be noticed that in Eq. ([8]) we neglect the contribution of the 
ion temperature to the sound velocity since typically Tj <^ Tg. However we 
keep the ion temperature in the nonlinear term ~ V^|E|^ since, as we will 
see in Sec. [3l it is this term which is responsible for the soliton stabilization. 
In the rhs of Eq. ([8]) we also neglect term ~ — ?7,o(a)V^|Ep/M, which is 
small compared to the contribution of the Miller force. Moreover Eq. ([8]) we 
take into account that (E^) = 6|Ei|^ while averaging over the time interval 



3. Cubic-quintic nonlinear Schrodinger equation 



In this section we sliall derive the nonlinear Schrodinger equation for 
the amplitude of the electric field. This equation will be considered in a 
particular case of a radial plasma oscillation. We shall numerically analyze 
the characteristics of the corresponding Langmuir solitons. In particular, 
their stability will be examined. 

Let us suggest that the density variation in Eq. ([8]) is slow, i.e. ff^nldt^ <^ 
c^g'W'^n. In this subsonic regime Eqs. ([7]) and ([8]) can be cast in a single NLSE, 
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which has both cubic and quintic nonlinear terms. NLSEs analogous to 
Eq. ([9]) were examined previously in connection to the studies of the light 



bullet propagation in crystals [U 



We shall examine axially or spherically symmetric plasma oscillations, 
i.e. E = Eeri where e^ is a unit vector in radial direction and E is a, scalar 
function. Introducing the following dimensionless variables: 
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we can represent Eq. (|9]) in the form, 
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which contains no dimensionless parameters. Here c? = 2, 3 is the dimension 
of space. 

One can check by the direct calculation that the following quantities: 
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(the plasmon number) and 
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(a) (b) 

Figure 1: The analysis of Eq. (fTTI) in 2D case, (a) The function N{X). (b) Examples of 
the numerical solitons "ipoi^) for different A. The solid line corresponds to A = 0.44, the 
dashed line to A = 0.42, and the dash-dotted line to A = 0.36. 



(the Hamiltonian) are the integrals of Eq. fill I) . Here ^^2 = 2tt and ^3 = Air 
are the sohd angles in two and three dimensions. 

We shall look for the solution of Eq. f lTT|) as ip{x,t) = e'^^'^ipo^x), where A 
is a real number meaning the dimensionless frequency shift. By the proper 
choice of the phase we can always make the function tpQ to be real. The 
corresponding ordinary differential equation for the function ipo is solved 
numerically using the MATLAB program. It requires an appropriate initial 
guess function ipg{x). The guess function is taken to have the Gaussian form 
il)g = Axexp(— a;^/2cr^), with the parameters A, a, and A chosen to minimize 
if ([ED at the constant A^ (^. 

Firstly, we analyze the stability of axially and spherically symmetric soli- 
tons by plotting A^(A) dependence. It is shown in Fig. [T]^a) for 2D case and 
in Fig. Mjo) in 3D case. One should notice that in 2D case dN/dX > in 
a quite broad range of A. Thus according to Vakhitov-Kolokolov criterion 
(VKC) [l2| this kind of 2D solitons is stable. In 3D case the accuracy 
of calculations is significantly lower that that in 2D situation. To build a 
smooth A^(A) curve in Fig. [2]^a) the least squares method was used since the 
points on this plot, obtained with numerical simulations in MATLAB, have 
rather big spread, especially at large A. Meanwhile one can see that unstable 
and stable solitons coexist in 3D case. In Fig. [2](a) we get that dN/dX < 
at A < 0.26. Applying VKC we conclude that this branch corresponds to 
unstable solitons. However, if A > 0.26, we obtain that dN/dX > 0, which 
corresponds to stable solitons. 




(a) (b) 

Figure 2: Panels (a) and (b) are the same as in Fig. [T] but correspond to 3D case. In 
panel (b) the solid line corresponds to A = 0.41, the dashed line to A = 0.25, and the 
dash-dotted line to A = 0.07. 



4. Application 

In this section we discuss tlie application of the results obtained in Sees. [2] 
and [3] for the theoretical description of a natural atmospheric plasmoid called 
a ball lightning (BL). 

BL is a glowing object appearing mainly during a thunderstorm. Despite 
the existence of numerous BL models (see, e.g. Ref. |l3|), it is likely to be 
plasma based phenomenon. The most spectacular BL property is its long life- 
time, up to several minutes, which is untypical for an unstructured plasma. 



In Refs. |14j . Il5| we put forward a hypothesis that BL is a radial oscillation 



of plasma and discussed both classical and quantum implementations of this 
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kind of oscillation. Note that analogous idea was considered in Ref. 

It is interesting to mention that spherically and axially (snake-like BL, see 



Ref. 17|) symmetric oscillations of charged particles do not have a magnetic 
field. Thus this system do not lose energy by electromagnetic radiation. 

The main problem of the existence of a plasma oscillation is its stability 
against the collapse jj]. Various mechanisms for the arrest of the collapse 
of a Langmuir wave in plasma, such as accounting for electron nonlinearities 
and additional quantum pressure of electrons, were studied p, Ig, l7| . 

In the present work we put forward another mechanism for preventing 
the Langmuir wave collapse based on the interaction of the ions EDM with 
the oscillating electric field of a plasmoid. We assume that ions in plasma 
are diatomic and nonpolar. Nevertheless it is allowed for an ion to acquire 
an induced EDM in the electric field. Thus this situation corresponds to a 
atmospheric plasma mainly composed of N^ (^ 78%) and O^ (~ 21%) ions. 



It is known that neither N^ nor O^ possess EDM because of the symmetry 
reasons. That is why spatial Langmuir sohtons studied in the present work 
are important for the description of stable atmospheric plasmoids. 

For the stability of a soliton its size should be greater than td - other- 
wise thermal electrons will cross the soliton leading to the development of 
a turbulence. For a rough estimate we shall take the typical ion temper- 
ature Ti = 300 K and suppose that the electron temperature in plasma is 
Te = 10^ K. One can see in Figs. [U^b) and El^b) that in dimensionless vari 



ables the typical soliton size can be up to 10. Using the results of Ref. [18 



for the polarizabilities of a neutral nitrogen molecule aij ~ 10~ cm~ and 
Eq. ( ITU]) , we get that for uq > 10^° cm~'^ a plasmoid becomes stable. This 
value is close to the density of an ideal gas at standard temperature and 
pressure, the Loschmidt's number, n; = 2.7 x 10^^ cm~'^. 

The analysis of the soliton stability is based on the perturbation theory 
which requires that, cf. Eq. fITOl) . 

, 15 Ti 1 
1287r2 Te (noAa)2 ^ ' 

Using the above chosen parameters of plasma and A ~ 0.4 in Figs. [U^a) 
and [2](a), we get that the condition flT4l) is satisfied if uq > 10^^ cm~^, 
which is one order of magnitude bigger than the previous constraint on uq. 
Nevertheless, one should notice that the nonlinear ponderomotive force fpoi 
in Eq. (jl]) will be defocusing in any case, providing the soliton stability. 
If condition ^\M is violated, one should just analyze the complete set of 
nonlinear hydrodynamic equations for plasma rather than use the simplified 
perturbative treatment. 

Finally can estimate the typical size of a plasmoid requiring that -Rgff > 
ro. For no = 10^° cm'^ and T^ = 10^ K, we get that i?cff > 10~^ cm. Such 



nano-sized plasma structures were discussed in our previous work [IJ]. Note 
the obtained tiny core of an atmospheric plasmoid can explain the fact that 
sometimes BL passes through small holes or cracks in dielectric materials 
without changing its shape and without destroying surrounding materials. 

5. Conclusion 

In the present work we have studied stable spatial plasma structures pos- 
sessing axial and spherical symmetry. The stability of soliton-like plasmoids 
against the collapse is provided by the defocusing interaction of the induced 
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EDM of ions with the rapidly oscillating electric field. Note that an ion was 
supposed to be diatomic without static EDM. In Sec. [H starting from the 
complete set of plasma hydrodynamic Eqs. ([1]), (|2]) and (jl]) we have derived 
the basic nonlinear equations for the envelope of the electric field ([7]) and 
for the perturbation of the ion density (|8]). In Sec. |3]we have reduced these 
equations to a single cubic-quintic NLSE fITT]) . written in the dimensionless 
variables f lTU]) . Then Eq. f lTT]) was solved numerically, cf. Figs, [^b) and[2](b). 
Applying VKC, we have found that mainly stable spatial solitons can exist 
in 2D case, Fig. [U^a), whereas in 3D case both stable and unstable solitons 
are present. Fig. |2]^a). 

In Sec. IHwe have discussed a possible application of the obtained results 
to the theoretical description of a long-lived atmospheric plasma structure, 
BL. It was found that in a realistic atmospheric plasma, mainly made of N^ 
and O^ ions, one can expect the existence of stable plasmoids described in 
the present work. Note that Nj and O^ are diatomic and can have only 
induced EDM, i.e. the our approach is valid. 

We have found that, for the preventing of a turbulence development, i.e. 
to guarantee the plasmoid stability, the density of background ions should be 
greater than no = 10^° cm~^, which is very close to the Loschmidt's number. 
The slightly greater value of uq compared to rii can be explained by the 
fact that we have used the polarizabilities of N^ as for a neutral nitrogen 
molecule, which in fact can differ significantly. Nowadays no measurements 
of N^ polarizabilities have been made. 

We have also obtained that the typical plasmoid size is small, Res ^ 
10"^ cm. Such nano-scaled plasma structures were discussed in our previous 



work [IJ] where we applied the quantum mechanical approach to explain 
their stability. 

In our estimates of the plasmoid characteristics in Sec. HI the temperature 
of the ion component of plasma has been chosen as 300 K, which is an 
appropriate value. However, to provide the reasonable plasmoid properties, 
the temperature of the electron gas had to be high, Tg = 10^ K. We remind 
that the electron temperature in a linear lightning discharge can be slightly 
higher than ~ 3 x 10^ K |19|]. Thus to generate this kind of plasma structures, 
described in frames of our model, one has to implement rather significant 
heating of electron gas. We can notice that a plasma with nn = 10^° cm~^ 



and Tg = 10 K may be created by a nano-second laser pulse |20|. This fact is 

important for the laboratory generation of spherically symmetric plasmoids. 

Finally we can also notice that, taking into account the small energy 



~ 10^ J of an object described in frames our model, one may consider this 
kind of plasmoids as proto-BL, i.e. an atmospheric plasma structure at its 
initial stages of evolution. Such a plasma oscillation can appear when strong 
rapidly oscillating electric field, probably associated with a linear lightning, 
is applied to a natural nano-sized spike, creating a plasma with a very high 
electron temperature. Under certain conditions, i.e. when other nonlinear 
effects become important, this proto-BL can then be transformed into a 
glowing object identified as BL. 

Appendix A. Polarization of a nonpolar diatomic gas in an exter- 
nal electric field 

In this Appendix we shall derive the permittivity of a gas of a nonpolar 
diatomic molecules in an external electric field. 

The collective response of plasma results in the concept of "dressed" par- 
ticles. It means that the Coulomb interaction of a charged particle is replaced 

by the Debye-Hiickel potential, ^exp [—-^)- If the plasma density is quite 
high, i.e. the Debye length is short, we may suppose that the electric charge 
of an ion is quite perfectly screened by surrounding electrons. In this case the 
Hamiltonian of a diatomic ion possessing an axial symmetry in an external 
electric field reads 

H=^ + U, (A.l) 

where M^ is the vector of the angular momentum of an ion perpendicu- 
lar to the ion axis, / is the moment of inertia of an ion, U = — (pE) = 
— E^ (q;_i_ sin^ 9 + a\\ cos^ 9^ is the potential energy of a polarized ion in an 
external electric field, and 9 is the angle between the ion axis and the electric 
field direction. 

On the classical level the thermodynamical properties of a gas can be 
calculated on the basis of the canonical partition function Z = z , where J\f 
is the total number of ions. The reduced partition function z which includes 
the rotational degrees of freedom of an ion has the form. 



d^M_L exp 



{27rhyj -—^-^y 2ITi 



dr^exp 



■j]2 

— (q;_l sin^ 9 + a\\ cos^ i 
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(A.2) 



where dQ = 27r sin 6d6 is the sohd angle differential. 

Calculating the integral over the angular momentum components, we can 
express z^t in the following form: 

Zrot = ^T, exp (frttxl z' , z' = j dxe«^' = 1 + ^ + |^ + . . . , (A.3) 

where ^ = E^Aa/Tj. The polarization of the gas in an external electric 
field can be calculated on the basis of the expression for the free energy 

F = —Ti In Z, as 

1 fdF\ ^ f , , 4 (AaE)2\ 

where V is the volume of a gas. Here we account for the decomposition of z' 
in Eq. COll . 



Basing on Eq. (]A.4p we obtain the permittivity of the ion component of 
plasma 

e = l + 87rn,(^(a) + A(^|p!^, (A.5) 

which was used in Sec. [2] to derive the ponderomotive force acting on ions. 

Finally we mention that our calculations are based on the assumption 
of the constant electric field whereas in Sees. [2] and [3] the electric field is 
supposed to oscillate with the high frequency ~ Wg. Using the results of 
Ref. [l8J one finds that the typical frequency associated with the polarizabil- 
ity of a molecule is ~ 10^^ Hz, which is several orders of magnitude greater 
than plasma frequencies in realistic plasmas. Thus the approximation of the 
constant electric field is valid. 
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